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Abstract 

Global existence for a system of nonlinear partial differential equa- 
tions (PDE) modeling an isotropic incompressible viscoelastic ma- 
terial is proved. The structure of the PDE is derived through con- 
stitutive assumptions on the material. Restriction on the size of 
the initial displacement and velocity for the model is specified in- 
dependent of the size of the viscosity of the material. The proof 
of global existence combines use of vector fields, local energy decay 
estimates, generalized Sobolev inequalities, and hyperbolic energy 
estimates. 
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1. Introduction 

This paper is concerned with the Cauchy problem for three dimensional isotropic 
viscoelasticity. The equations for the motion of the viscoelastic material are de- 
rived from constitutive assumptions. The resulting system of equations is a hybrid 
between the equations for a viscous Newtonian fluid and a Cauchy elastic material. 
These equations can be viewed as a quasilinear symmetric hyperbolic system with 
a parabolic perturbation of a certain form arising from the viscosity. Studying the 
system in this way, this article presents a result concerning the global stability of 
the motion of this material. 

The main theorem states that the system has global solutions in time given 
sufflciently small initial displacement and velocity, and furthermore, the size of 
these initial data are independent of the magnitude of the viscosity. The methods 
used for this result include generalized hyperbolic energy estimates, generalized 
Sobolev inequalities, and local energy decay. The parabolic term is avoided for the 
most part and docs not provide additional decay during the course of the proof. 

Although many treatments of elasticity use Lagrangian coordinates, the equa- 
tions will be studied as a first order system in Eulcrian coordinates with constraints 
as in [14j. Using these coordinates, the form of the constraints is particularly useful 
and the portion of the linear operator due to the pressure waves is easily seen to be 
higher order. Written in Eulcrian variables, the resulting system is a mixed type 
PDE, an interpolation between the Navier-Stokes equations and incompressible 
isotropic elastodynamics. 

The proofs in this article utilize previous techniques and insight from the study of 
3D quasilinear wave equations and elasticity. Global existence cannot be expected 
from a general symmetric hyperbolic system. For certain types of "genuine nonlin- 
earity" John proved blow up for the quasilinear wave equation in which the second 
derivatives of the solution blow up in even for arbitrarily small initial data [4] . 
In turn, he showed a blow up result for elasticity using similar strategies on equa- 
tions containing these types of "genuine" nonlinear terms 5 . As a complement to 
John's "genuine nonlinearity" , Christodoulou 3 and Klainerman [J concurrently 
identified a "null condition" for the quadratic terms of the wave equation. When 
all the quadratic terms satisfy this null condition, small data global existence can 
be proved. Klainerman's techniques involve proving strong dispersive estimates 
and using of Lorentzian and Gallilean invariant vector fields in the energy without 
directly estimating the fundamental solution. 

Analogous global stability results for nonlinear elasticity could not immediately 
be realized because the equations contain multiple wave speeds in the linear oper- 
ator causes a lack of Lorentz invariancc. For clastic materials which are Galilean 
invariant and isotropic, works by John and later Klainerman and Sideris [8] 
showed "almost global" existence for small initial using a smaller set of vector 
fields. In (llj Sideris went further, introducing a null condition (similar to that of 
Klainerman and Christodoulou) for both the shear and pressure waves present in 
compressible elasticity. He showed that the shear waves inherently satisfy the null 
condition while the pressure waves generally may not. Using these null conditions 
and the previous vector field methods he proved global existence for compressible 
elastic materials close to equilibrium without estimating the fundamental solution. 
Independently, Agemi achieved the same result using null conditions and estimation 



of the fundamental solution J^. Blow up for elasticity with large initial data was 
shown by Tahvilar-Zadeh in |15j . 

Small data global existence for incompressible elasticity was proven by Sideris 
and Thomases using the incompressible limit in |12| and directly in [14j . Both 
results use vector field methods and dispersive estimates to prove global existence. 
In the case of isotropic elasticity only the inherent null condition for shear waves 
is needed because the pressure waves do not appear, so no additional assumptions 
were made on the material in either result. In 13J, Sideris and Thomases obtained 
dispersive estimates in the framework of a more general result on local energy decay 
for symmetric hyperbolic systems. Although systems with certain dissipative terms 
are handled in these estimates, the techniques cannot be directly applied to the 
viscoelastic case because the viscosity does not appear in all components of the 
PDE. 

Previous results on the Cauchy problem for incompressible three dimensional 
viscoelasticity by Lei, Liu, and Zhou f9], Lin, Liu, and Zhang [ lOJ , and also Chen and 
Zhang i2j use the parabolic structure of the equations to achieve global existence. In 
each article, the authors study the Oldroyd B model for viscoelasticity, equations 
with a structure nearly identical to the equations studied herein. Their global 
existence proofs involve writing the PDE in terms of a "special quantity" which is 
a specific linear combination of the deformation gradient and the Eulerian gradient 
of the velocity. This change of variables reveals additional dissipative structure 
which is then used to prove small data global existence. For these results, the 
magnitude of the initial data must be small relative to the viscosity. Consequently, 
convergence to the equations for elastodynamics cannot be studied in this way 
because the initial data vanishes with the viscosity. 

To study global solutions for viscoelasticity from the perspective of 3D quasi- 
linear hyperbolic PDE, the paper is organized as follows. The PDE are derived 
from constitutive assumptions on the material. After introducing some necessary 
notation and stating the main theorem, local energy estimates similar to those of 
[T5] are computed. Projections are then defined on the eigenspaces of the symbol of 
the hyperbolic operator, and a null condition for shear waves is stated in terms of 
the projections. As usual the pressure is shown to be bounded by nonlinear terms, 
and weighted Sobolev inequalities are stated. A new addition to these estimates is 
combined with a Hardy-type inequality to give the key inequality (|8.12|) . Following 
bootstrapping of the local energy estimates and computation of the energy identity, 
all of the pieces are pulled together to prove global existence. 

2. Derivation of Equations 

First, assume our substance is a continuous distribution of matter at rest and 
that it fills three dimensional space so that each point X — {Xi^X2tX^) G R'^ 
corresponds to a point in the substance. These points X are called material or 
Lagrangian coordinates. If we deform our material in a differentiable manner to 
a different configuration at a later time t, we call these new coordinates spatial 
or Eulerian coordinates and label them x{t,X). The map x : R x R-^ Y{? 
is an orientation preserving diffeomorphism. The inverse map X{t,x) returns the 
material point that has deformed to spatial coordinate x at time t. The deformation 
gradient is the matrix 



4 



with inverse 

^ dx 

Both H and F must have positive determinant because of our assumption that x 
preserves orientation. We denote 

D= — = {DuD2,Ds) and V = — = (^i, ^2, ^3) 

as well as Dt for Lagrangian time derivatives and dt for Eulerian time derivatives. 
Using the chain rule, we have the relations 

= F^dp and Dt = dt + v ■ V 

where v = Dtx. Here and in what follows, we use the Einstein summation notation 
where summing over repeated indices is understood. 

Constitutive assumptions must be made to further specify material properties. 
For an elastic solid the equations of motion are typically derived from a variational 
problem, however the dissipative motion of a viscoelastic material prevents us from 
obtaining our equations in this way. Therefore, we begin with the force balance 
laws expressing the conservation of mass and momentum: 

(2.1) Ap + V-fp = 

pDtv -V-T = pf. 

where p{t, x) is the density, f{t, x) is the external force (which we assume to be 0), 
and T(t, x) is the Cauchy stress. If we write these equations with Eulerian time 
derivatives , we have 



(2.2) dtp + v-Wp + V ■vp = Q 

pdtv + u • Vu - V • T = 0. 

We first assume the material is incompressible. This ensures the density is 
constant and transforms the conservation of mass equation into the constraint 

(2.3) V • t; = 0. 

This assumption also means volumes are conserved under the motion so that 

(2.4) detF = deti? = 1. 

For simplicity, we will assume the density is unity for the rest of the paper. 

The Cauchy stress T encodes the internal self-interacting forces of the material. 
For a viscoelastic material T depends on the pressure p{t,x) (linearly), the defor- 
mation gradient F and its Lagrangian time derivative DtF. By the chain rule, we 
can write the ODE 

DtF = WvF 
so the Cauchy stress can be written as 

T=-pI + f{F, Vv) 

where / is the identity matrix. Assuming objectivity (Galiean invariance) of our 
material implies that T depends on Vf through the rate of strain tensor D'j = 
^{djV^ + diV^). If we assume this dependence, which determines the internal fric- 
tional forces within the material, is linear the Cauchy stress takes the form 

(2.5) T = pI + uoD + f{F). 
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Here the tensor T contains the information about the clastic properties of the 
materiaL If we assume the elastic forces come from an isotropic and objective strain 
energy function W{F) as in [Ullllj, then we have 

W{FU) = W{F) = W{UF) 

for all proper orthogonal matrices U . The isotropy and objectivity assumptions 
imply that W depends on F through the principal invariants of the strain matrix 
F'^ F. The Piola-Kirchoff stress, S, is the Lagrangian version of the Cauchy stress 
and is defined as 

dW 

In this case, since our material is incompressible and p = 1, T is related to W by 

■ dW 

(2.6) ^pW;^^^ = ^p^p'^^^ = 

We will assume the material is stress free at the identity, i.e. 

f (/) = S{I) = 0. 

This restriction rules out the Oldroyd B model. 

Taking Eulerian divergence of (|2.6p and using djF^ = (a consequence of in- 
compressibility) , we have 

flW B^W 

(2.7) [V . f{F)Y = Wg^iF)) = d^i^d^Fj: ^ At{F)Fld,Fj: 
where A = ^ \s the elasticity tensor. A is symmetric, i.e. = A^^, because 
(2.8) 



dF 

dFjOFf ~ dFfdFj 



We remark that this form of the elasticity tensor is consistent with the notation in 

mm- 

We will not make any further assumptions about the nonlinear terms of these 
elasticity tensors, however, to be consistent with linear elasticity theory, we impose 
the Legendre-Hadamard ellipticity condition 

(2.9) Ai^(/)e,ac^V'>0 

for all ^, a; € S"^ , the unit sphere in R"^. This is a standard assumption on the form 
of the linear operator ruling out the incompressible Navier-Stokes equations. Under 
the objectivity and isotropy assumptions this condition is equivalent to 

where 5 denotes the unit tensor and the parameters ci and C2 are positive constants 
representing the propagation speeds of the pressure and shear waves respectively. 
Because of the convenience of the relation 

we will write our PDE in terms oi H = F~^ defined above. Using the fact that 
FH = / we have 

d,Fl^ = -F-F^dpHl. 
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So we can write 

(2.11) A^^^{F)F^d,F^ = AL^(F)5^Ff9pF," 

= A%{F)FlHrFp^Fj: 
= -A%{F)FlF^F-FlHrd,H:^. 
If we add the null Lagrangian = c^iS^S"^ — can define 

(2.12) ia(i7) = A%{H-')[H-']l[H-^Y^[H-^UH-% + 

so that A{T) is positive definite. This does not change our equations since 

via the constraint (|2.10|) . A satisfies 

APJ^iH) ^ Arr^iH) 

and, also, for \H\ = \H — I\ sufficiently small, 

(2.13) M:^,{H)H^HI>cI\H\\ 
by continuity of A. 

We are now in a position to write our PDE in terms of H and v. First, by (|2.5p . 
(|27f|) and (l2lT|) we have 

(V • Ty = -d.psi + "^djidjv' + d,v') + (V • fy 

= -d,p + v^v' - A';j{H)HfdiHl 

where v = ^ and we used the incompressibility constraint V • u = 0. Our conser- 
vation of momentum equation becomes 

dtV + v-Vv + Vp + N^i^H) = vAv 

where N'"{Hy = A^j^{H)HfdiW^. 

Using the chain rule and the definitions of H and v we have the transport equa- 
tion 

dtH + vVH + H\/v = 0. 

Thus, the equations we have derived for a 3D isotropic incompressible viscoelastic 
material are 

(2.14a) dtH + vVH + HVv = 

(2.14b) dtv + v -Vv + Vp + N^iH) = i^Av. 

As we have already seen, the constraints we have accompanying these PDE are 

(2.15a) V-u = 

(2.15b) detiJ^l 

(2.15c) d,Hl=dkH] 

At this point we remark that the Oldroyd B model for viscoelasticity corresponds 
to the case where above 

47(i?) = {H-'y;{H-%{H-X{H-X- 
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Although this form of A cannot be derived via our method, the form of the equa- 
tions is the same and the constraints are precisely the same as they arise solely 
from incompressibility. The Oldroyd B system is typically written in terms of the 
deformation gradient F as 

(2.16a) dtF + v -WF - FVv = 

(2.16b) dtV + v-Vv + Vp-V ■{FF^) = iyAv. 

To verify equations (|2.14aP - (j2.14bp are equivalent to (|2.16a|) - (|2.16bp in this case, 
one can check using F = H^^ that 

We also note that the null condition which is discussed in section 6 is also satisfied 
in the Oldroyd B case. Therefore, Theorem 14.11 holds for the Oldroyd B model as 
well. 

Turning again to the general case, we linearize the system (|2.14ap - (l2.14bp using 
the notation H = H — I and v ~ v obtaining the equations 

(2.17a) dtH + Vv ^ N"{H,v) 

(2.17b) dtv + V -TH + Vp-iyAi^ N''{H,v) 

where 

{N^TjiH, v) = -v^'dpH} + Hppv' 

{Ny{H,v) = ~iPd,v^ - A'^{H)HfdiHl [A{H) A{I)y^m^„ 
(V • THY = A^\mHL = d,[clHl + {cj - cl){tr H)Sl]. 
Our constraints then become 
(2.18a) V • w = 

(2.18b) tr ij = i((tr - tr (H^)) + det H 

(2.18c) djHl ^ dkH] 

Notice that the incompressibility of the system implies that tr H is higher order in 
H. Because of this, wc will move the tr H portion of the linearity onto the right 
side of (j2.17bp . For simplicity, we will assume herein that C2 = 1. We caution the 
reader that our dot notation iJ, v does not denote derivatives, consistent with the 
convention in [13j . After a bit of rearranging we have the equations in the form we 
will work with most frequently: 

(2.19a) dtH + Vv ^ N"{H,v) 

(2.19b) dtV + V -H - vAv = N'"{H, v) - {cf - 1)M"{H) - Vp 

where 

{Ny{H,v) = -i^d.e - {Alj{H)HfdiHi, + [A{H) A{I)]^diHl,) 
(M^'yiH) = -^.[^((tr ij)2 - tr (H^)) + det H] 
with constraints (|2.18ap - (|2.18cp . 
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We also take this opportunity to write the system in the notation of [Mj with 

(2.20) LLf = dttl - AkdkU - (0, t/Aw) = N{U) - (0, {c[ - 1)M"{H) + Vp) 
where N{U) = {N" [U), N"" {(])). 

3. Notation and Preliminaries 

Throughout the paper we use the following notation without mention. Most of 
our norms and inner products will be in and most integrals will be taken over 
R'^, so we write 

MI = II-IU^(R3), 

•) = •>L2(R3), 

and 




The notation 

(c> = (c^ + l)^ 

is used as a substitute for the real number c when being bounded away from zero 
is important. 

We use the usual derivative vector fields 

d = {dt, 81,82,83) and V = {81,82,83) 
as well as the scaling operator 

S^t8t+r8r, 

its time independent analogue 

= r8r, 

and rotational derivatives 

n = (.T A V). 

Here, as usual, 8r = f ■'V. 

Because we are working with vector valued functions we will also need to use 

(3.1) n^{u) = {a,H + [v^'\H],n,v + v^'''>v), 

where 

1/(1) = 63 ® 62 - 62 (8 63, 

y(^' = ei «) 63 - 63 ® ei, 

and 

= 62 «) 61 - 61 (gi 62. 

The notation U = {H, v) and U = {H, v) is used frequently, where H, H are three 
by three real- valued matrix functions and v, v are three dimensional real valued 
vector functions. The arguments of these functions are nearly always suppressed. 
Occasionally we will write 

(3.2) QH ^nH+ [V, H] for e ® 
and 

(3.3) nv = nv + Vv for w e R^ 
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as notation for the components of flU. For scalar functions / we define 

(3.4) nf = nf. 

For convenience we use the notation 

T ^ (V,0). 

An exponent a on T denotes an ordered fc-tuple, (ai,...,Q!fe) £ Z'', where 
makes sense for 1 < i < fc. We denote 

and 

|a| — k. 

By TU we mean TiU for some i. 

The following identities and commutation properties of the vector fields are used 
as well. One extremely useful identity which we will use often is 

(3.5) V = u}dr--{uj 

r 

for Lo = ^. The commutator of any two Ts is a T. We note that T is a Lie algebra 
with the Lie bracket given by the commutator operation 

[T„T,] = T,T, -T,T,. 

We denote 

(3.6) A{'sJ)U ^ AkdkU ^ {S/v,SJ ■ H) 
and recall the definition (|2.20p 

(3.7) LU ^dtU + A{V)U +{Q,-vAv). 

Because of our assumptions of objectivity and isotropy, L commutes with T, that 
is 

(3.8) T"Lf/ = LT°'U. 

However, S does not commute with L. For a function / and n = 1, 2, 3, ... we have 
the identities 

(3.9) = 9(5' - 1)'7 
and 

(3.10) S-^A/^ ('"V-l)"^^(5- 1)^/. 
So, if we use the notation (5* — 1) = S* we have 



(3.11) s'^Lu ^[dt^ Aiyy^s'^u ^vA^i M-iT^^sw. 

One of our main concerns will be dealing with this commutator. 
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4. Spaces, Norms and Statement of Main Theorems 
Writing So = rdr, the space for the initial conditions is 
Hl'^ = {U^ {H,v) : R3 ^ (R3(g)R3)xR3 | S^T'^U e L^{R.^),\a\+a < a,a< 9}. 
with norm 

^ |Q|+a<<T 

Solutions will be constructed in the space 

ffp''^ = {U = (H, v) : [0, oo) X R3 ^ (R^ ® R^) x R^l 

gajo^jj g C'([o,oo);i2(R3)) for a<9,a+\a\< a}. 

Define corresponding norm to be 

||c/||^.,« = sup{(t)-^ E Ws^r'^uw} 

^i^^ |Q| + a<<7 

where 5 < 1 as in Theorem 14. II below. Given U G H^'^ define the energy by 

(4.1) E„^g[U]^^ J2 I (ip7W5'''T"ijfS'"T"ij/„ + IS-^T^wpj . 

\c.\+a<a "'R-^ 

As long as \U\ is small, as will always be assumed, we have 

(4.2) E,,e[U]^ E WS'^^Uf 

Q| + a<<7 

a<e 

via (|2.13p and the standard Sobolev Embedding Theorem in three dimensions. This 
equivalence will be used repeatedly in the sequel without further comment. The 
time independent analogue of E is 

(4.3) £aAU] = \ E / (i^T(ff)5o"T"i/f5o"T"ij4 + |5o"T"t;|2). 

Using these definitions, we state the main result. 

Theorem 4.1. Let Xq he an orientation and volume preserving diffeomorphism on 
R'^, and let vo be a divergence free vector field on R"^. Define 

Uo^{Ho,vo)^{VXo,vo) 

Uo - [HoM = {Hn-I,vo). 

Suppose that Uq G H'^'^^'^ with k > 17 and that 

£k,k-3[Uo] < C, £^-4.^^-4Uq] < e 

for uniform constants C and e. 

If £ is sufficiently small, then the initial value problem for I2.14a^ - I2.14h\l with 
initial data Uq = U{0) has a unique solution U{t) which satisfies the constraints 
h2.15a]) - \2.15c\) and is a member of H^''^^ . Furthermore, the magnitudes ofS^.K-siUo] 
and £i^-4^i^-4[Uo] do not depend on the size of the parameter v and the solution sat- 
isfies the estimates 



11 




a+|Q|<K— 4 



(4.5) £;„,«_3[C/(i)] + ^ / llV^^T^lpdr <C'£;,,„_3[C/o](i)' 




for all t G (0, oo), where C is a uniform constant and S < I. 

The proof of this result, which will be given in section 10, uses local energy decay 
estimates, Sobolev inequalities and hyperbolic energy methods. 

We remark that these estimates reduce to those of 3D incompressible elastody- 
namics (see [M]) as ly 0. 

5. Local Energy Decay 

In order to take full advantage of energy estimates, we need to derive inequalities 
that establish decay away from characteristic cones for solutions of the system 
(|2.19ap - (j2.19bp with constraints (|2.18ap - (|2.18cp . The isotropy of the system implies 
all of the necessary hypotheses for local energy decay which are discussed at length 
in [14] . The subtle new feature of our equations is appearance of viscosity term 
only in equation (j2.17bp . When = or if we artificially add a viscosity term, 
—i/AH, to (|2.17aP the following result reduces to [14]. 

Before we begin the proving the estimates, we must specify the regions in R x R-^ 
where the estimates will hold. These regions are dictated by the light cone = t 
associated with the linear problem. Define a smooth cutoff ^ : R ^ R which 
satisfies 



C = 1 on [0, 1 



and 



C = on [2, oo). 



Let m > 4 be a fixed integer. Define 



(5.1) 




and 




Theorem 5.1. Suppose U = {H,v) is a solution of 



dtH + Vv 



f 



(5.3) 



where f and g are forcing terms and that U also satisfies the constraints 12. 18a]) 
and i2.18c]) . Then the following bounds hold for n = 1,2,3... and exponent a — 
{ai,...,ak) whenever the right hand sides make sense: 
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(5.4) ^{||77(<)V^''T"{7|| + Wrji^tAS'^T^vW} 

n 

^ X!^''"^'^"^^^"^" + ll'S'°+^T"{7|| + ||VS'°T"[7|| 

+ ||S'"T"C/|| + ||r;iS'"T"/|| + \\T]tS''T°'g\\ + iy\\7]tW ■ S"'T"/||}. 



(5.5) Eill^t'^^'- ~ i^(V)]^"T"C/|| + \\-/,.tAS^T"v\\} 

n 

< ^{||fi^°T"ij|| + \\S''+^T°'U\\ + ||VS'°T°w|l 

a=0 

+ ||S'°T"w|| + ||7iS'"T"/|| + WjiS^T^gW}. 

Furthermore, the suppressed constants do not depend on the viscosity v. 

Proof. If we take derivatives 5'°T" in equations (|5.3p . use commutation properties 
of the vector fields (|3.8p and (|3.1ip . we get 

(5.6a) dtH + Vv^ f 

(5.6b) dtv + W ■ H - v/\v = g 

where {H, v) = U = S'^T^'U, f = S^T^f and 

(5.7) g = S-^T^s + ) (-l)""^'5'^T"t). 

Multiplying equations (|5.6ap - (j5.6b[) by t, moving the tdt terms to the right and 
rewriting them using the scaling operator gives 

(5.8a) tVv = rdrH - SH + tf 

(5.8b) tV ■ H - vt/\i = rdri> - Sv + tg. 

Equations (|5.8a|) - (|5.8bp will be our starting point for the derivation of the estimates 
(113) and lEM- 

Interior Estimate, (r < — ) 

To begin we multiply equations (|5.8a|) - (j5.8bp by rj, apply || • |p to both equations 
and use the triangle inequality to obtain 

(5.9) WvtVif < 2\\vrdrHf + MlvSHf + Mlvtff 

(5.10) \\TjtV ■ Hf + \\7]iytAvf - 2(77tV • H, r^vtAv) 

< 2\\'nrdrif + iW-qSiW^ + AW-qtgW^. 
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where above and in the future we use the notation (•, •) (•, ■)l^{-r3). The main 
calculation is estimating the cross term in (|5.10p . Taking divergence of (|5.8a|) . using 
commutation properties, and multiplying by rjiy we get 

(5.11) T]i'tAi = r^v{rdrV ■ H -V ■ SH + tV ■ f). 
Thus, 

(5.12) -2{r]tV ■ H,r]i'tAi) 

= -2{tjW ■ H, r]v{rdrV ■ H ~ V ■ SH + tV ■ f)) 
= -2{rjtV ■ H, r^vrdrV ■ H) + 2{r^tV ■ H, rjvV ■ SH) 

~ 2{7jW ■ H,r]iy\7 ■ f) 
^ {II) + {12) + (/3). 

For any g = 1,2, 3..., 

(72) > -gi^^\\7]V ■ SHf - -\\r]tV ■ 

Q 



and 



(/3) > ^gy'^WijiV ■ fW" --Wi^tV ■ HW" 



Q 

by Young's Inequality. 

Estimating {II) is a more significant calculation. For all p = 1, 2, 3... 
(71) = -2vt J Tj^dkHlxjdjdpH; 
= -2vt I ij'x,d,{^\S/ ■ H\') 

>-2^tJr^\r,'\^\V-H\' 

>-Avt j 7/max|?7'||V • 77p 

= -2 j {2vTnax\ri'\\V ■ II\){rit\V ■ II\) 

> "-llTytV • Hf - 4:giy^ma.x\r]'f\\\/ ■ Hf. 
Q 

Here we used integration by parts, ^ < 2 on supp 77, and Young's Inequality. The 
value of g will be determined later. Adding estimates, we have 

(5.13) -2(?7<V ■ 77, f]vtl\i) = (71) + (72) + (73) 

> — ||77tV • 77||2 - 4£<i/2^ax|77'|^||V • 77f 
Q 

- gv'^Wi]^ ■ SHf - gv^W^qtV ■ /f . 
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Combining ()5.10p and (I5.13P gives 

(5.14) (1 - -)||?7tV ■ Hf + Ihi/iAwf 

9 



< 2\\T]rdrv\\'^ + 4||77S'^f + 4||?7<gf + igi^'^ max \rj'\ \\V ■ i/f 
Now we compute 

(5.15) WrjtV -Hf^ J if^dkHldpHl 

= -J 2W<'^^^9p#;- J v't'HidkdpH; 

>-2j rjtmW\\H\\V-H\+ J 2rj,^'t''^HldkH; 

>-2m J \r]'\\H\{T]t{\\/H\ + \V ■ H\)) + J T]h^\VH\^ 
> -Wmax|?7'p||7?f - ^{\\r]tVHf + \\T]tV ■ H\\^) 

where we used integration by parts twice, constraint (|2.18cp . and Young's Inequality. 
Estimate (|5.15p implies 

(5.16) -\\TjtVH\\^ - 9m^max\f]'\^\\H\\^ < \\r]t\/ ■ H\\^. 
5 

If we combine (|5.9p . (|5.14p and (|5.16p . we have 

(5.17) {l-^)%tVUf + \\rj,.tAif 

g 5 

< 2\\rjrdrUf + A\\TjSUf + A\\ritgf + ^r]tff + 9m^ msiic\Tj'f\\Hf 
+ Aqv"^ max |r;'|^||V • Hf + gi^^hV • SHf + Qu'^htV ■ ff. 



Using W-qrdrfjW^ < Ih^^^f/p < :^\\ritVU\\'^, and taking = 8, we can absorb 
2\\r]rdrUf on the left: 

(5.18) (i_^)||^(t)v;7f + h^.fA^f 

< 4||,?5t/f + A\\i^t~gf + AW~ff + 8i.2||^iV ■ /f + hVt/f 

+ W max|7y'n|7?f + 32t/2 max |?7'|^||V • ijf + 8i^2||77V • S'i^lp 
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which we can write without exphcit constants as 

(5.19) ||7?(t)VS^"T"{7|| + Wr^iytAS^T^iW 

< ullVS^+^T^HW + WS^+^T^UW + \\WS''T°'U\\ + ||S'°T"{7|| 
+ ||?/tS'"T"/|| + ||77tS'"T"5|| + i^\\T]tV ■ S-^T"/!! 

+ \\rj,ytAj2(''){-'^r~'Sn"v\\ 

3=0 ^-^^ 



after recalHng the definitions of U, H, v, /, and g. We pause here to derive the 
corresponding exterior estimate. 

Exterior Estimate, (r > -^) 

Rearranging (|5.8ap - (|5.8bp . multiplying by 7, applying || • |p to both equations, 
using the triangle inequality and adding the resulting estimates together yields 



(5.20) hrdrU - 7<A(V){/||2 + W^ti'Aif - 2(7<V • H - -frdri, jiytAi) 
<2\\jSUf + 2\\jtff + 2\\^tgf 

Again, most of the work is calculating the cross terms. The simpler of the two is 



(5.21) 2(-/rdrV,-fi'tAv) ^2vt / YrdrV^'Av 



-Aut I jj'^x,d,v'dkv' - 2vt I ^''dkffdkff 



2vt j 7 XjdjdkV^dkV^ 



> -Aut I j\j'\'j^\drvf - 2ut J j^\Vv\^ 



^2vt / ^^'^xjdkv'dki' + ivt / -i^\VvY 



> -Sut I 7|7'||a^f}p - Avt I 7l7'|-^|Vu|^ 



mr^ 



+ vt I 72|Vvp. 



> -12i^tmax|7'| Jj\Vv\'^, 



where we used integration by parts, and ^ > 2 on supp 7'. We pause to integrate 
by parts 



16 



(5.22) - / 7|V{)|^ = - / ^dkv'dkv 



{t)r J 

/Tfl f 
l7'l^l«l|VS|-y 7l«l|Aw| 

> -m(t)~i max I7' I j \v\\\7v\- j \v\\-i/\v\. 
Now (|5.21l) . (|5.22p and Young's inequality imply 

(5.23) 2{-irdrV,-ivt/\v) > -12i/i max |7'| J ■^\Vvf 

> -12t/max|7'|^ / | Vi)| - 12 max |7'| / \v\\-fiytAv\ 



> - max |7'|2 f 72e||«f + VSf ) - ^hi/tASf , 

for any positive integer g. 

To estimate the other cross term in (15.201) . we start by recalling (|3.5p 

V = Ujdr - - {uj AQ) 

r 

to get 

(5.24) \I-H=[udr--{LoAn)]-H 

r 

^uj-drH - -{u An) H. 
r 

Using (|5.8ap we have 

(5.25) tiv ■ drH = tio ■ (-Vu + -SH - -/) 

r r r 

t ~ ~ 

— — Lu ■ Vv H — uj ■ SH uj ■ f. 

r r r 

Combining (|5.24p . and (|5.25p gives us 



(5.26) -2{jW ■ H,-fiytAi)) 

t 

= -2(7— w • Vw, ji'tAv) - 2(7-tj • SH, -fvtAv) 
r r 

~ 1 . 

+ 2(7— w • /, -ivtAv) + 2{-ft-{uj An) ■ H, jiytAv) 
r r 

= {El) + {E2) + {E3) + {E'i). 
Using Young's inequality and the fact that — < 1 on supp 7, we can bound 



{E2) + {E3) + {E4) below by 
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(5.27) --W^vthif - Qm^W^SHr + Wltff + hnHf) 

Q 

Now we only have one term left to bound: 



(5.28) {El) = -2i^t^y 7V^2a;fe9fcu'Au* 



So far, we have only used integration by parts. Continuing our estimate we repeat- 
edly apply the facts ^ < 1 on supp 7 and \drv\ < \yv\ to get 

(5.29) (£^1) > -Gvt^m^l^ J -/\j'\{ty\mry^\Vv\'^ + J j'^{mr)'^\Vv\'^^ 
-2vclt^m^ I j^{t)-^{mry^\Wv\^ 



> -6i^im^(^max|7'| J 7|Vw|^ + J 7^|Vu|^) - 2vm^ j -i'^\S/v\'^ 

> -22vtm^ui&yi\"f'\ j -f\Vv\^ 

> -22TO^i/max 17'p y |{j| | Vu| - 22m^ max |7'| j \v\\-ivt/\v\ 

> -max|7'|2(242m'*^?||^}f + i^^||V{)f ) - -\\-ivt/\if 

Q Q 

for any positive integer g. In the last two lines we used (j5.22p and Young's inequality. 
Combining ((03| . (jO?)) . and ((5?28)) with g = 10, we have 

(5.30) -2(7iV • H - •jrdri, •jiytAi) 

( jyTTL^ ^ 

> - max |7'|^(720 + 2420m'*)||wf - max \-f'\^- ^|| VuH^ 

5 

- lom^hsnr + mr + hmn ^ii7^<A«f . 

Inserting this into (|5.20p gives 
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(5.31) hrdrU ~ "ftA{V)Uf + i||7ti/Awf 

< 2||7t9f + (10m2 + 2)(\\^SUf + ||7</|P + h^Hf) 

+ max |7f (720 + 2A2Qm%vf + max |y pi^^ || V{;f 

5 

which can be written without expUcit constants as 

(5.32) W-irdrS^T'^U ~ 7tyl(V)S"'T"C/|| + \\-itv ISS^T"" v\\ 

< ||7tS'°T"5|| + ||^°+iT"C/|| + ||7tS'°T"/|| + ||i7S'"T"ij|| 

g-i / \ 

+ ||5°T"f/|| + ||V^°T"i;|| + hi^tAj2 I ) (-l)""-''^'T"i>|l . 

Induction on a. 

Estimates (|5.19p and (|5.32p are analogous, corresponding to (|5.4p and (|5.5p 
spectively. For each, we must deal with the final term on the right hand side, 
will use an inductive argument to handle these terms. 

To prove (|5.4p . first assume a = 0. Using (|5.19p . we have 

(5.33) ||77(t)VT"f7|| + ||77z/<AT"t)|| 

< \\VST°'H\\ + ||S'T"{7|| + |lVT"f7|| + \\T"U\\ + \\TitT"f\\ 
+ \\rjtT"g\\ + htV-T^f\\. 

For a = 1 

(5.34) ||77(t)VS'T"{7|| + \\riiytAST"v\\ 

< \\VS^T°'H\\ + ||S'2T"f7|| + ||VS'T"[/|| + ||S'T"J7|| + \\T]tST°'f\\ 
+ \\7]tST"g\\ + \\rjtS7 ■ ^T"/|| + Hryi/t AT"i)|| 

which is bounded by 
1 

(5.35) ^{||VS"^+iT"H|l + \\S''+^T"'U\\ + ||VS'"T"C/|| + |1^"T";7|| 

+ ||77i5°T"/|| + ||77i5°T"g|| + ||7/tV • 5'^T"/||} 

after applying (|5.33p to the final term of (I5.34p . 

More generally if we assume (j5.4p holds for n ^ k, then we have 

(5.36) ||,7^iA^('"V-l)'^-^5^T%|| 

k 

< \\7ji^tAS^T"'v\\ 

j=o 

k 

< ^{||VS'^+iT"ij|| + IIS-^+iT^C/ll + ||VS'^T"f7|| 
j=o 

+ WS^T^ilW + \\j]tS^T°'f\\ + WTjtS^T^gW + \\TjtV ■ S'^T"/||} 
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which impUes (using ()5.19|) ) 



(5.37) 



]{t)VS''+^T"if\\ + \\riiytAS''+^T"'v\ 



k+l 



+ ||S'"T"[/|| + ||r;iS"'T"/|| + \\r]tS''T°'g\\ + W-qtV ■ S'°T°/||} 

Combining this estimate with our inductive hypothesis gives (15.41) for general n. 
(|5.5p foUows from an identical argument using (|5.32p . 

□ 



6. The Null Condition and Spectral Projections 

For u! G S'^ let A+ = 1,A_ = —1 and Aq = be the eigenvalues of A{uj) = 
AkOJk with Ak defined as in (|2.20p . Denote the orthogonal projections onto the 
corresponding eigenspaces by F+{uj),¥^{uj) and Po{uj). Using the formula 



(6.1) 

we compute 
(6.2) 



^+{uj)U = -{{v + Huj) (E)u;,v + Hu) 



(6.3) 

and 
(6.4) 



^_(w)[/== -{-{v ~ Hoj) ® u},v - Hoj) 



M^)U = {H(I -uj(E,uj),0). 



For notational simplicity we will suppress the argument lu in our calculations. 
Using these projections we utilize a calculation that appeared in [Hj . 

(6.5) \{X,t~r)P,d.jS''T"U\ 

= \V,{tA{uj) -rI)djS''T"'U\ 
< \{tA{uj)~rI)djS''T°'U\ 
= \{tAk - ruj''I)uj^djS''T'^U\ 

= \{tAk - rLo''I)iu^dk + ir!fc,)^''T"J7| 

r 



1 



uj\tA{V) - rdr) + {tAk - ruj''I)-nkj 



< |(tA(V) -ra^)S'°T"f7| +C 



^ + 1 
r 
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Now we rewrite estimate (|5.5p as 

3 n 



(6-6) J2 E E^ll^^^'^ - r)P,a,^"T"C/|| + \\-/,ytAS''T"v\\} 

j—l i a— 

n 

< ^{llf^^^T^f/ll + ||^''+iT"{7|| + ||V^°T"i;|| 

+ ||^'^T"{>|| + ||7t5''T"/|| + ||7i5'^T"g||} 

where lo — 'm each expression above. 

The null condition is an inherent property of the quadratic coefficients in the 
nonlinear terms of our PDE. Formally, it states that the quadratic terms are linearly 
degenerate causing them to decay enough to prove global existence. Without a null 
condition, we may have a genuine nonlinearity (see [5]) meaning a global existence 
proof would not be possible via hyperbolic energy estimates. The way we use the 
condition can be seen in the energy estimates (|10.18p - (|10.24|) . 

In our work with the null condition it is convenient to define 

Pi{ljj)u ~ U — {u,Lu)ji3Ll! 

for u G R'^. Pi(w) is the orthogonal projection onto spanjijj}-'- and P2(ct') is the 
orthogonal projection onto spanjo;}. Again, the argument lu will be implied in 
calculations. 

The elasticity tensor A satisfies a null condition restricting the quadratic inter- 
action of shear waves. More precisely, for -Bp™ — gj^i , A''™- 5^ we have 

(6.7) i?;,7,"(/)c.zc.„c.„efi)e^2)C(3) = 

for LU e S*^, ^(i) G span {w}^, i — 1,2, 3. Using our projections we can write 

(6.8) B^7fe"(/)^zc.™c.„(Pi^(i))P(Pie(2))^'(PiC(3))' = 

for all LU E S'^ , ^(i) G R'^. The tensor A^^S"^ clearly satisfies the null condition since 

(6.9) 47(/)<5I?c.zc.,„c.„(Pie(i))P(Pie(2))^'(Pie(3))' 

(6.10) = i;,™(/)u;,c.,„(Pie(i))^(Pie(2))^(c^, (Pie(3)))R3 
= 0. 

If we write in terms of the elasticity tensor A we have 

If we differentiate term by term, the most difficult term comes when the derivative 
falls on A 

(6.12) -^[A'^dAFmF^Fl^F::! ^ ^^Ali,{F)^\F^\FlF]Fl^ 



-±^AP^l{F)F^F-FlF]FtF^^ 
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which is equal to 

d 

when evaluated at the identity. The proof that this tensor satisfies (|6.7I) can be 

p 



found in [TT] . When the derivative falls on Fp we have 



d 

which equals 



(6.13) - Ali{F)^[Fl]F]FtF:^ = AI-Uf)FIF^F]fI^ F^^ 



at the identity. This satisfies (|6.7p by (|6.9p . A similar argument holds for differen- 
tiation of the other F terms. 

For bookkeeping purposes we define the weighted LP' norms 

(6.14) X,,e{U)=Y.Yl E ll(A,<-r)P,a,5-T";7|l 

(6.15) 2.,«(C>)-EE E ll7(A.t-r)P,,9,^''T"f/|| 

t 7 |<:.|+a<,T-l 

and 

(6.16) ^.fi{V)^ E h(i)V^''T"{7|| 



| + a<<T-l 
a<8 



These quantities will be used to link our local energy decay estimates to the energy 
estimates via our bootstrapping lemma to follow. 

7. Bound for the Pressure 

To eliminate the pressure term in local energy decay, we will bound it by the 
nonlinear terms and iV" using the equations and constraints for Fl. 

Lemma 7.1. Let a = 0, 1, 2, 3..., and a he a suitable exponent for T. Then for p 
satisfying equations h2.19a\) - i2.19h\] and constraints V2.18a]) - li2.18c\] . 

IIV^^T^pII < ||5"T"iV''([7)|| + ||5°T"M^(ij)||. 

Proof. If we apply derivatives 5"'T" to equation (|2.19ap - (|2.19b|) and rearrange we 
can write 

(7.1) VS-^T^ = S-^T^Vp 

= S"'T"(iV''(C/) + (1 - c\)M"{H) -V ■H + vMi - dtv). 
We proceed by taking the divergence of the equation to obtain 

(7.2) AS"'T> =V • S'°T"iV''([/) + (1 - c?)V • S"T°'M"{H) 

- V ■ (^°T"V ■H) + uy ■ S^T^Av ~ V • S^T^dtv. 
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Using the constraints (I2.18a|) - (|2.18c|) we can simplify the last three terms: 

(7.3) -V • (5''T"V ■ H) ^ -d^d^S^T^H] 

= -d^S^'T'^djiv H 

and 

(7.4) t/V • S'°T"Ai-; - V ■ 5°T"9tu = {lyS^T^A - 5°T"c»t)(V • v) = 0. 
So (|7.2p simplifies to 

(7.5) AS-^T^p =V • S'°T"7V''(C/) - c?V • S''T'^M"{H). 

Using the notation = 5°T"(7V^)*({/) - c?S'°T"(Af^)^(ij), we can estimate 

(7.6) ||VS'"T>||2 = -(A5''T>, S'"T» 

= -(V-$,5''T» 
= ($,VS"'T» 
< ||<i>|| . ||V^«T>1| 

which implies the result. □ 

8. SoBOLEV Inequalities 

We need several Sobolev inequalities beyond the standard embedding H'^iTi?) ^ 
i°°(R3) to do our bootstrapping and energy estimates. 

Theorem 8.1. For A, e R, A G [0,2], / e C^(R3), r = \x\, and p = \y\, we have 
the following: 

(8.1) r\f{x)\<Y, \\dp^V{v)\\%y\>r) X E \\^yiy)\\LHlyl>r) 

\a\<l \a\<2 

(8.2) r{X,t-rmfix)\<Y, 1^ - p}d,nViy)\\iil^^^^^ 

\a\<l 



\a\<2 



(8.3) 1/(^)1 <E \\P'%^''f(y)\\LHlyl>r) 

\a\<l 

X Y: ||/-^f^"/(y)|li(^(|,|>.). 

|q|<2 

Proof. Estimates (|8.ip and (|8.2p are proven in ^llj Lemma 3.3. 
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To derive ()8.3p first assume g{x) G C(f (R'^). Then using Cauchy-Schwarz and 
calculus facts we estimate 

(8.4) / \g{rcj)\*du; < r [ \dpg{pLu)\ ■ \gipL,)\^dpdiu 



< 



p-'\dpg{y)\-\giy)\'dy 

p'>r 

1/2 / ^ \ 1/2 



\Jp>r / \Jp>r 



Using Gagliardo-Nirenberg, we have 

(8.5) ||5(MlUa(s.)< J2 m"giru)\\%.)MrLo)\\%%2^ 

\a\<l 

SO 



(8.6) / p'('-'^\g{y)fdy= / p'^-'+'ng{pu;)fdu;dp 

Jp>r Jr JS^ 

< 



|a|<l 

Putting ([Q]) together with ([8^ . we get 



I |i^"g(pc^)pdc^ I ]dp 

^sup / \g{puj)\'^duj 

p>r Js^ 

X E / p'^'-'^\n'^g{y)\'dy. 



1/2 



(8.7) (^^^ |5(rc^)|^dc.) < p-''\d,g{y)\'dy^ 

Y: I p'^'-'^n-g{y)?dy 

Combining (|8.7p with the isopcrimetric Sobolev inequality 

(8.8) |/(x)| < l|f^"/MllLMs^), 

|a|<l 

we complete the proof of (|8.3p . □ 



Lemma 8.1. ("^^ard?/; J'or / e C5"(R^), r = |a;|, and p = |?/|, 
(8.9) \\p''f{y)\\ < \\dpfiy)\\. 
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Proof. The result is implied by the estimate 

„-l Y'^..M|2 / „-2| /•/..m2 



\f{puj)\'^dujdp 

'0 JS^ 

ndp{p)\f{pi^)fdujdp 

Jo Js^ 



=2 y p-Vp/(y)l • \f{y)\dy 

<2\\d,f{y)\\-\\p-'f{y)\\. 



□ 



Corollary 8.1. Let ui — r — \x\, and U — {H,v) satisfy constraints i2.18a\) - 

mm- If 

Ea+ \a\-\-2.,a \a\-\-3,a {U),'^a+\a\+3,a{U) < OO, 



then 



(8.10) (r)|S"'T"[/| < ^ ||S"'T"+'3f/|| 

l/3|<2 

(8.11) {r)(X,t~r)^/'^\¥,S''T"U\<Y^ ||S'"T"+^f/|| 

\p\<2 

+ ^a+\a\+3,aiU) 

(8.12) ||(f)5'^T"[/|U^(,>o) < J2 WS'^r'^^^UW 

l/3|<2 

+ ^a+|c.|+3,a(f^) 

(8.13) r^/^i^ . (5«T"ij^)| < ^ ||5"T"+'3[7|| 

|/3|<2 

+ ^ ||r5"T"+'5M^(ij)|| 

I/3|<1 

(8.14) r3/2|^ . (^^x"^)! < J2 ll^"T"+^f/||, 

l/3|<2 



where rj is defined as in I15.1\) . E as in l{4.1\ l, X and ^! as in ^6A^ and 116. 16\) . and 
M" as m lilM] . 

Proof. Inequalities (|8.10p and (|8.1ip are proven in [TT] Proposition 3.3. To prove 
(|8.12p notice that on supp 77 we have 

(8.15) \{t)S'^T^U\ < |77(t)5"T"{7| + ||(i)^''T"{7|Uoo(o<,<i). 
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Applying (18. 3p to the first term on the right with A = 1 and / = ri{t)S°'T°'U gives 

(8.16) \v{t)S''T"U\<J2 ||p"^apfi^(ry(t)S"^T°C/)||i/2 

I/3|<1 

X J2 \\p-'^^v{t)S''T'^U)\\'/'. 

I/3|<2 

Applying (|8.9p to both terms on the right and using Young's inequality we have 

(8.17) \r,{t)S''i:"U\ < J2 \\dl^^iv{t)S''T''U)\\ 

I/3|<1 

+ J2 \\d,n^ms''T"u)\\. 

I/3|<2 

Using the fact that flrj ~ rjil, we can bound the first term on the right 

(8.18) Hni'ms^T^mw 

I/5|<1 

< E {\\dpWs''r''+^u)\\ + ||ap(77(i)ap5'^T"+^t7)||) 

I/3|<1 

^ i\\{tr'v"S''T''+'^U\\ + l|ry'V5"T"+^[/l| 

I/3|<1 

+ \\r]{t)V^S^T''+'^U\\) 

< Y ||^'^T"+'^f/||+*a+|a|+3,a(f/). 
I/3|<2 

The second term on the right of (|8.17[) is bounded similarly. To complete the proof 
of (|8.12p we use {t) < (r) on {0 < r/ < 1} and (|8.10p on the second term on the 
right in (|8.15p to obtain 

(8.19) mS^T'^Uh^^o<n<i) < \\{r)S^T'^U\\L^ < ^ \\S''T'^+^U\\. 

\I3\<2 

Proving (|8.13p and (|8.14p begins with applying ([Q]) with A = i and / = r^/'^uj-g 
to get 

(8.20) r^/^lc ■g\<Y Wp-'^'d^n^p'/'c ■ 5)llii(|,|>,.) 

|q|<1 

X Y ||p-/^f^"(p3/2..5)||^/4„>.) 

|a|<2 

<E llp-^/^5p(p3/2c..f]"5)ll 

|q|<1 
|a|<2 

<E iip^^5pfivii+ E 

|q|<1 |q|<2 
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Using the identity Ij3.5|) we can estimate the first term on the right 
(8.21) \\roj'drn°'g'\\ = ^ ||r(9, + r-^{n A uj)i)n"g'\\ 

\a\<l \a\<l 

|a|<l |a|<2 

In hght of (11201), taking g = S^T^v yields ^^J^ after applying (|2.18ap while taking 
g = S'^T'^Hu) and applying (|2.18bp and (l^TScl) gives 

(8.22) 

E ||r-f2^5,(S'"T"ij]wJ)|| < E {\\^^^S''T°'H\\ + \\rVt^d,S''T°'H]Lo^\) 
I/3|<1 l/3|<i 

< E (||Vfi''S'"T"ij|| + ||rf7'5^"T"A/^(ij)||), 
l/3|<i 

which implies (|8.13p . □ 
We remark that the constraints (|2.18ap - (|2.18cp are only necessary to prove (|8.13p 

and (mm). 

9. Bootstrapping and the Energy Identity 

Before we can begin bootstrapping we need a technical lemma to deal with terms 
where multiple derivatives fall on the elasticity tensor A. A similar lemma appeared 

in m]. 

Lemma 9.1. Suppose U E H^'^ with cr > 3. Set a' = [(t/2] + 2. Suppose 
Ea-'.a'[U] < 1 and \U\ < S for all t G R+, with S sufficiently small. If d is any 
positive integer then for f : (IV^ (g) R-^) x R-* R'' satisfying \f{U)\ = 0{\U\p) at 
the origin we have the pointwise estimate 

\S^T'^f{U{t,x))\< E \S'''Tl^'U{t,x)\- ■ ■\S^^i:'^''U{t,x)\ 

bi+...+bp<b 
/3ll + --- + l3pl<l|3 

forb+\(3\ <cr, |/?| <e. 

Proof. Using the chain rule, write 

(9.1) S'''T^[/(f/)] = E E f^^^''\U)S'''T'^'U ■ ■ ■ S^^+'T^^^+'-U 

j<b bi+..^+b-^^=b 

where /^"^ denotes the nth derivative of / with respect to U. At most one derivative 
may exceed order [it/2] because + b < a. By the standard Sobolev embedding, 
commutation properties and (|4.2p . we have 

|S"=T'^f7| < \\W^S''T'U\\ < ||(5 + 2)"V2T'-"f/|| < E„,^„,[U] < 1 
whenever c+ < [cr/2]. By the mean value theorem, 

\f(i+'^\U)\<\U\P-=-\ j+k<p 
for \U\ < 1. Our result now follows from (|9.ip . □ 
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In order to make use of our local energy decay estimates we use our assumptions 
of the smallness of the initial data and the hyperbolic nature of the system to bound 
our weighted norms So-,e,^cr.e, and Xa-,e defined in (|6.14p - (|6.16p by the energy. We 
accomplish this by bootstrapping the nonlinearity. 



Lemma 9.2. Fix k > 15 and ii = k — 4. Suppose k — 1>(7>1, fi>9>l and 
U{t,x) solves equations \2.19a\) - i2.19h\l and satisfies constraints !i2.18a]) - \2.18c\) . 
As long as e]i{'^'\U] is sufficiently small 



(9.2) ^^AU)<eI'^^^,^,[U] 

(9.3) '^^fi{U)<El'^^,^^[U] 

(9.4) X,AU)<eI'^^,^,[U] 



each hold whenever the quantities on the right make sense. 



Proof. Combining Theorem l5.1[ Lemma FTT] and equations (|2.19ap - (|2.19bp gives us 



(9.5) S.,,([/)=^^ {\\l{Kt-r)l?,d,S'^T'^U\\ 

t i |Q| + a<<T-l 

< {||j7^''T"{7|| + ||S'°+iT"C/|| + ||VS"'T"i;|| 

|q| + „<„_1 

+ ||S'°T"w|| + \\-itS"'T'^N"{H,v)\\ 

+ \\^tS''T'^[N%H,v) - {cl - l)Af«(ij) - Vp)||} 

|c,| + a<a- 
a<9 + l 

+ Y HII'S"'T"iV({7)|| + ||S'"T"i\/^(ij)||} 

l£»| + a<<T-l 
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and 

(9.6) ^.AU)= E \\vmS^^"U\\ 

Q| + a<<T-l 
a<0 

< J2 {||V5"'+iT"ij|| + ||S"'+iT"C/|| + ||V^°T"f/|| 

l€,] + a<c,-l 

+ ||S'°T"C/|| + Wr^tS^T^N^iH,!})]] 

+ ||77<5°T"(iV''(ij,i;) - (cj - \)M"{H) - Vp)\\ 

+ \\T^tV ■ S''T°'N"{H,v)\\} 

< J2 \\S''T''if\\+ J2 Wvty ■ S''T°'N"{H,v)\\ 

|c.|+a<<T+l |Q| + a<<T-l 

a<e+i a<e 

+ t{\\S''T"NiU)\\ + \\S''T"M"{H)\\} 

q| + <x<<7-1 

a<e 

1 /2 

Our strategy will be first to bound by Ej^^ for low energy levels (cr < 
^ — 1), then use this result to bound E^^e by -E^+i g^i for low energy levels. We 

^ 1/2 

will use these results to bound 5'^,^ and ^(j^e by -Eg-^i g_^_i at high energy levels 
(/^ < cr < K — 1). 

Assuming cr < /it — 1 we begin by estimating using Lemma |9.1[ splitting into our 
interior and exterior cutoffs and applying (|8.10p to obtain 

(9.7) t\\S''T°'N{U)\\ 

< Y {||77t|S"=T^f7| • |VS"'T''{7||| 

o|+„<^_l c+d<a 
a<8 k| + |f|<|c. 

+ ||7t|S"=T'-"f/| • |VS"*r''f7|||} 
^ E E {h(^)I^^T^t/|-|V5'^T^{/||| 

Q| + a<„_l c+d<a 
a<8 k| + |iS|<|c 

+ ||(r)|S"=T'^[/| • \VS'^T^U\\\} 

Q| + a<„_l c+d<a 
a<8 k| + |0|<|c 

+ ||(r)S"=T'^[/||Loo||V5''T''{7||} 

< Y Ell^[U]{\\ri{t)VS''T°'U\\ + \\VS''T°'U\\}. 

„i+„<„_i 

1 /2 

As long as i^^;^ [?7] ^ 1 is small enough, this implies 

(9.8) i|l^°T"7V(f/)l| < CeI[^[U] + e*,,e(f>) 

|Q|+a<„_l 
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where e <iC 1. By a similar argument we have 

(9.9) i|l^"T"Af«(ij)|| < CEyj[U] + e*,,e(f>). 



-a<cr-l 



The remaining term from (j9.6p is 

(9.10) J2 lk/iV-5'"T"7V^(ij,w)l| 



|qI+„<„_1 
a<fl 



-HZ! {||?7(i)|V5"T^[/| • |V5'^T''[7||| 



|q| + o<<t-1 c + d = a 



+ ||??(t)|S"=T'^L/| ■ |V^5''T''C/|||}. 

If we temporarily assume WLOG that c + |<^| < d + then since c + |<;| + 3 < 
f + 3 < /i the first sum on the right can be estimated as follows 

(9.11) Yl ll'7(^)|V5=T^{7| • |V^'*T^{7||| 

^ Z Z llV5"T'^t/|U-h(i)V5'^T*{7|l 



3l + a<o--l c+d=a 



c«| + a<<T-l c + d = a 
Q| + a<<7-l 



< 



e«'.,9(f/) 



where again e ^ 1. We will need to work harder to bound the second sum in (|9.10p . 
The exceptional case is when = a+ |a| = cr — 1 because there are too many 

derivatives to take the weight ri{t) with the term. Instead we use (|8.12l) to get 
the estimate 

(9.12) ||?7(t)|;7| ■ |V'5"T"[/||| 

<h(t)f7|U-||V25"T"f/|| 

<( E l|T^t^||+*3,o(f/)) WV^S'^T^iJW 

\l/3|<2 / 

When we are not in the exceptional case we have d+\-d\ <ct — 2so that 

(9.13) ||?7(i}|S"=T^C/| • IV^S'^'T'';/!!! 

<\\S''T''ij\\LA\ri{t)V'^S'^T^U\\ 
<e«'a.e(f/). 

This leaves us with the bound 

(9.14) Y ll^iV-5«T"7V^(ij»|| <Ci?y+',_,[[/]+evl,,,e(f/) 

\a,\ + o.<a--L 

a,<a 
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Altogether, dM]), jM]), , and ((9111) give the estimate for ct < ^ - 1. 

Keeping the assumption cr < ^ — 1 we have from ()9.3p . (19. 5p . (19. 8|) . and (|9.9|) 

(9.15) S,,e(f/)< 5] ||5"T"[/|| 

Q| + a<<T 

a<e+i 

+ t{|iS'°T"iV(C/)|| + ||S'°T"M^(ij)|l} 

Q| + a<<T-l 

<Eli\.e+i{U] 

proving (|9.2p for cr < /i — 1 . 

For the rest of the proof we will assume ^ < a < k — 1. As in (|9.7p we have 

(9.16) ^ i||S'"T"iV({7)|| 

c|+a<<T-l 

Q|+a<cr-l c + d<a 
a<8 |5| + |tf|<|a 

+ \\jt\S''T'U\ ■ \VS'^T'^U\\\} 

^ E E {Ms^^'ui-ivs^r^uiw 

\a\+a<cr-l c + d<a 
a<e |5| + |0|<|q 

+ ||7i|S'"T^{/| • \VS'^T^U\\\}. 
For terms where < ^ — 3 we have 

(9.17) ||?7t|5"T^[7| • |VS"*T''[7||| + ||7i|S'"T^[/| • |V5'^T''{7||1 

<||S'"T^;7||Loo||7yiV5'^T''{/|| + ||(r)V5''r''f7||Loo||S'^T^;7|| 

where we used \I'^_i^0(C/) < E^{^(U) < 1 in the final inequality. On the other 
hand, for terms where d + > /i — 2, we know c + j";! < 4 allowing us to get the 
following bound: 

(9.18) ||77t|S"=T'^{7| • \VS'^T'^U\\\ + ||7i|S"=T'^J7| • \VS'^T'^U\\\ 

<||S"=T'^{7||Lo=||77iV5'^T''{7|| + || (r)S"=T'-"f/||Loo || V5'*T''{7|| 

<eI'][u]{^,Au) + eI'M) 

Combining (|9.16p . (|9.17p and (|9.18p we have the estimate 

(9.19) A\S''T'^N{U)\\ < CE'J^,^,^,[U] + e^>,AU)- 

\a\ + a<a-~l 
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And by a similar argument 

(9.20) t\\S^r-M"iH)\\<E'J^,^,^,[U]+e^.AU). 

a| + a<<T-l 

a<e 

To finish the proof of (j9.3p we need to estimate the remaining nonhnear term in 
()9.6p . Starting with (|9.10p we bound the first sum by assuming WLOG that c+|<^| < 
(i+ |t?| because both terms contain a gradient. Then since c+|?| + 3< ■| + 3<^we 
can follow estimate (|9.1ip . As for the second sum on the right in (|9.10p we again use 
(|8.12p combined with (|9.3p in the exceptional case when d + li)] = a + \a\ = a — 1: 

(9.21) ||7/(t)|;7| ■ |V25'^T"[/||| 

<||7y(i)[/||Lo.||V2^°T"Lf|| 

<| J2 IIT'^f^ll +*3,o(f>) I ||V2S"^T"C/|| 

\l/3|<2 / 

< {e1^o'[u] + eI!,'[u])eI/^\,[u] 

When not in the exceptional case we have d+li?! < a— 2 which means if c+|<^| < fj,—2 
then 

(9.22) |l77(t)|S'^T'^{7| • \V^S'^T^U\\\ 

<\\S''T'U\\L^\\'n{t)V^S'^T^U\\ 

and if c + |<;| > ^ — 1 then d + j-i?! < 4 giving us 

(9.23) ||77(t}|S"=T^C/| • iv^s'^'T'';/!!! 

<||5"T'^f7||Loo||ry(t)v25''T''{/|| 

<E'J+i,eP]^<^AU) 

after again applying (|9.3p for low energy. 

Taking ([CT|) . (PT^ - (p:^ together we derive estimate (H^. In light of 

(|9.15p we have also proven (|9.2p . The bound l|9.4p follows immediately since 

□ 

Before beginning our energy estimates we state the energy identity which can be 
obtained by a standard calculation for quasilinear symmetric hyperbolic systems. 
Notice we have not commuted the vector fields with the laplacian in the viscosity 
term. We will leave this commutation until the end of the energy estimates where 
we do induction on 9. 
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(9.24) dtE^4U]-iy ^ S''T°'v'S''T°'Av' 

a+|a|<,^ J 
a<0 

+ J2 J2 I Cb,a,0^o.(s''T°'v'S^Tf^vPdpS''T''v' 

+ 0{ Y H I l'5"T"[/| • \S''TPu\ ■ \S'^T^U\ ■ |V5"T^f7|)) . 

Here we have separated the quadratic terms which will be handled via the null 
condition and constraints from the higher order terms. We also note that we have 
used Lemma inn] and the smallness of the low energy to get the higher order terms. 



10. Energy Estimates 

Now we prove the main theorem by estimating the time derivative of the energy 
dtEcr^e\U] for high (cr = k) and low {a = k — 4 = fi) energy levels. Using these 
estimates we will do a finite induction proof on 9 to complete the result. We do 
not make any assumptions about until the inductive portion of the proof. 

Proof of Theorem \4-l\ 
High Energy 

Beginning by assuming a = k we estimate very roughly using (|9.24p . Sobolev 
embedding, and crashing through with absolute value to obtain 

(10.1) dtE^,g[U] - 1^ J2 f S''T"v'S''T°'Av' 

a<e 

^ H Y [ l'^''T"{/| • \S''T'^U\ ■ |V5"T^{/|. 



33 



We split our integral using the cutoff functions 77 and 7 

(10.2) J |S"'T"[7| • |S'''T''f7| ■ IVS^T'tll 

= J 7|5'^T"[/| • IS'^T'^f/l • IVS^T'U] 

+ J ri\S''T°'U\ ■ \S''T'^U\ ■ IVS^T'UI. 

To estimate the exterior integral we consider two cases and use the fact that 
r> on supp 7. If \P\+b< |^| + c, then |/3| + 6 < [f ] < ^ - 2, so 

(10.3) j 7|^''T"[/| • \S''T'^U\ ■ |V5"T^{7| 

< J TOr(i}-i|S'"T"{7| • |S'''T'3{/| • |VS'=T^[/| 
(t)-^ J {r)\S''T''U\ ■ \S''T^U\ ■ |V5=T^f/| 



< 

< (i)-i(r)5''T^{/||Loo / |5''T"f7| • |VS"=T^{7| 



Here we used the weighted Sobolev estimate (|8.10p . If + 6 > |^| + c, then 
kl+c< [^] + l <yU — 2, sobya similar argument we have the same bound in this 
case. Thus, 

(10-4) J2 J2 J2 [ llS^^^U] ■ \S''T^U\ ■ |VS'"T'={/| 

On the interior region we also consider a few cases and this time we use (|9.3p 
to recover a factor of (t)^^. The exceptional case is when |<j|+c= |a| + a— 1, 
and + 6 = 1 where we need to use (I8.12p in tandem with (|9.3p . In this case, we 
estimate 

(10.5) J 7?|S"'T"C/| • IS^T/^Ul ■ |V5'=T'^{7| 

< J (i)-i|5°T"[/| • |77(t)5^T^{7| • |V^"T^C/| 
<(<)-! II (t)5''T/'{7|Uo.(,>o) J l^'^T^f/l • |V5=T^;7| 
J2 \\S'T^+'U\\ + ^t+m+3AU))E^AU] 



< „ 

< 



<{ty'El/^,'[u]E^Au] 

<{t)-'Ei/mE.Au]- 



A second possibility is |/3|+fo < |<^|+c < |a| + a — 2 in which case |/3|+6 < [^] < ^ — 2 
giving the estimate 

(10.6) J r;|5°T"f/| • |S'^T^[/| • \VS''T'U\ 

< J (<>-i|^"T"C/| • \S''i:^U\ ■ \r,{t)VS'T'U\ 
<{t)-^\\S^Tf^U\\L^\\S''T"U\\ ■ \\rj{t}VS''T'U\\ 

using (19. 3|) and the standard Sobolev embedding. FinaUy, if |/3| + 6 > |<^| + c, then 
kl + c < [f ] - 1 < At - 5, so 

(10.7) J 77|5'^T"f7| • \S^T'^U\ ■ |VS"=T'^[/| 

^ J {ty^\S''T°'IJ\ ■ \S''T>^IJ\ ■ |?7(t)V5"T^{/| 
<{t)-^\\r]{t)yS^T'U\\L^\\S^T"U\\ ■ \\S''T^U\\ 

<{t)-'El('Au]E,Au] 

after using (|8.12p and (|9.3p . Taking estimates (|10.5p - (|10.7p together, we have 

(10.8) J2 J2 J2 [ vlS^T^il] ■ \S''T'^U\ ■ \VS''T'U\ 

<{t)-'El/'AU]E.,e+i[U]. 
Altogether our high energy estimates, (|10.4p and (jlO.Sp . give 

(10.9) dtE^AU]-!^ /'5°T"^;^S'°T"A^;^ 

^ J2 J2 J2 [ \S'"^'^U\ ■ \S''T^U\ ■ \VS''T'U\ 

a+|c|<» b + c = a c+k\=ia+\a\ 

<{t)-'El/'AU]E.,o+i[U]. 
In order to complete the induction on 9 later on, we will also need the estimate 

(10.10) dtE^^i,+i[U] - V /"^°T"t;'S"'T"Az;* 

a+|c|<K J 

<{t)-'El/l[U]E^,,+,[U]. 
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Considering (|10.4p with — fi + 1, we only need to show 

(10.11) H f r]\S''T"U\-\S''T^IJ\-\VS''T'U\ 

In terms where a < /z, we may apply estimate (|10.8p to achieve this bound, so we 
only consider the case when b + c = a = fi+l. Terms where b and c are each at 
most /z are estimated using (jlO.Sp - p0.7p . We are left with two exceptional cases, 
6 = /i + 1 and c = /i + 1. If 6 = /i + 1, then c — and < 3, so using (|8.12p and 
(jO)) we get 

(10.12) J 77|^°T"{7| • \S^T'^U\ ■ \VS''T'=U\ 

= J 77|5^+iT"[/| • \S''+^T'^U\ ■ |VT^{7| 
<(t}-i||5^+iT"[/|| • WS'^+'^T'^UW ■ ||r/(i)VT^{/||Loo 
<(t)-ii?«,^+i[C/](4(f [[/] + ^rAU)) 

A similar argument can be used when c = /i + 1. This shows (jlO.lip holds, so the 
estimate ()10.10p is valid as well. 

Low Energy 

For the low energy {a — fi) we need to rely on the precise structure of (|9.24p 
combined with the null condition (|6.8p to obtain decay on the exterior 

region. We split each integral on the RHS of (|9.24p into two integrals using our 
cutoffs rj and 7. On the interior region, we will easily obtain the desired decay by 
estimating roughly as in (|10.ip 

(10.13) J r;|^°T"C/| ■ IS-^T^f/l • |VS"=T'^C/| 

<(i)-2y" (<)27^|S'"T"{7| • \S''T'^U\ ■ \VS'=T'U\ 

J l'5"'T"{7| • |(t)S'^T'';7| ■ \r]{t)VS'T'U\ 

<{t)-^\\{t)S'T^Uh^^^>o)\\S''T'^U\\ ■ hmS'T'^UW 

<{t)-'{ E ll^''Tf''+'''^ll+*''+l/3|+3.b(f>))<|„|,jC/]*c+|.|+i,c(f>) 

\-/}\<2 

Here we used (f8J2l) . ([Q]) . and the fact that either 6 + < [f ] or c + |^| < [f ]. 
The exterior estimates begin with the computed energy from (|9.24p : 
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(10.14) J2 {[l[DnAl'^]iI){ld^v'^S-r^Hi^S''T^Hf 

^ _i_ I ^, I ^' ^ " 



i+|q|< 



J2 E Cb,a,p,c.S''T^H^lS^T^vPdiS''T'Hl,} 



et'^T^ c+|<,-|#a+|Q| 



J 7i^7(iI)^"T"ij^,a,ijfS"^T"z;* 

-L^ — r. _ ( I _ I / . ( I _ I 



.^J_w— ^ „1 1^1 / _ 1 \ ^ \ ^ 



We start by estimating the final term using r > — on supp -r. llS.lOp and the fact 
that either 6 + |/3| < [f ] or d + \d\ < [f ] 



(10.15) y 7|S"'T"[/| • |S'^T^f7| • |S"'T''[/| ■ |V^'=T'^f7| 



<{t)-^\\{r)S''T>^U\\L^\\{r)S''T^U\\L^\\S''T''U\\\\VS^T'U\\ 
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The convective terms from ()10.14p are estimated using ()8.14|) in tandem with (13. 5|) 
For example 

(10.16) J -/AlJ{H)S''T''Hi,S''T'^v'd^S''T'Hf 

7|5°T"{7| • \uj ■ {S''T'^v)\ ■ IS/S^T'UI 



+ {ty^ J 7|(r)5'^T"C/| • \S^T^U\ ■ \nS''T'-U\ 
+ {t)-^\\{r)S''T°'U\\L^\\S''T''^IJ\\ ■ WnS^T^UW 
<{t)-'/'E,,o[U]E'J][U]. 



The other convective terms 



and 



J jS''T°'v'S''Tf^vPdpS''T'v' 

have the same bound using a similar argument. 
We are left with the cubic terms 

(10.17) J 7[i:'f;i^™](/){ia„w''^''T"ij,45''T"ijf + diH^;S''T°'H^J''T°'v' 
E E / 7a,a,/3,«i^7(/){5"T"ffi^5^T^ijfa,^=T^ 

-L.-— ^ _ 1 I _l / . 1 I _ I 



which can only be estimated using the null condition. Each individual term above 
can be written in the form 



for some quadratic form Q where each of a+ la], b+ |/3|, and c+ |<^| + 1 are no larger 
than /X and ((•, •)) = (•, •)r30r^3xr3. If we decompose using the spectral projections 
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dOll-lEll) and (ISH), we have 

(10.18) J j{{Q\diS''T''U,S''T^U),S''T°'U)) 

K,L,Me{+-fl} 

= j i{{Q\^+diS''T''U,^+S^T^U)J+S''T°'U)) 

+ Yu [ l{{Q\^KdiS''T''U,f'LS^Tf^U),f'MS°-T°'U)) 

= j 7(((3'(P+c^i<9r5'"T^[/,P+S'^T'^f7),P+S"'T"f7)) 

- j 7((Q'(P+r-i(cj A n^S^T^il ,¥+S^T^U),¥+S''T°'U)) 

+ Y [ l{{Q\fKdiS''T''UjLS^TPu),¥MS''T°'U)) 

(K,LM)^( + , + , + )•' 

(10.19) ee(N1) + (N2) + (N3) 

We bound (N2) by using {t) < (r) on supp 7 and (|8.10p 

(10.20) J 7r^i|f]S"=T^{/| • \S''T^U\ ■ |S'"T"{/| 

< J -f{r){t)-^\nS''T'^U\ ■ \S''T^U\ ■ |5''T"f/| 

<{t)-^\\(r)nS''T'U\\L^\\S''T^U\\ ■ ||S"^T"f7|| 
<{tr^E,^g[U]E^,e[U] 

since c + |<^| + 3 < k. 

To estimate (N3), assume WLOG K +. Then because 

l<r{t)-' <{XKt-r){t)-' 

on supp 7, we have 

(10.21) J j{{Q\FKdiS''T'U,FLS^T'^U),FMS''T°'U)) 

< J 7(t)-^/2^r)(AKt-r)i/2|P;^V5"T^{/| • |S'^T'^{/| • \S''T°'IJ\ 
<(t)-3/2||(r)(AKt-r)i/2p^VS"=T^[/||Loo||S'^T'3[/|| • ||5''T";7|| 

after using ([8lT|) . . and c + |<;| + 5 < k. 
The remaining term, (Nl), is 
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Since each term in ()10.17|) has a different Q, we will carefully estimate one of these 
terms and leave out the details for the remaining four terms. The sum from the 
second line of (|10.17p has a quadratic form which gives the formula 

(10.22) J j{{Q'{F+ujidrS''T'-U,¥+S''Tf^U),F+S''T°'U)) 

= J j[D^^A';^]iI)u;i{S'T^'v + S''T''Hlj)''uj„ 

X [(Pi + P2)(S"'T"{. + S''T°'Huj)]P 

X [(Pi + P2){drS''TH + drS''T''Huj)y 

after decomposing using the formula (|6.2p and using our projections Pi and P2. 
This product yields eight terms, one of which is 

J 7[i5^it7](/V/[Pi(5''T^« + S'Tf'Hu;)]'' 

X [Pi(5''T"w + S'"T"ijw)]Pw„ 
X [PiidrS'^TH' + drS^T^Huj)]^ 
= 

by condition (|6.8p . The remaining seven terms all have at least one factor of P2, so 
each is bounded by 

where a+|a|, c+li^l < /i are generic exponents. Here we used the fact that the 

number of derivatives on each term in (|10.22p is less than n and [9^/1 < | V/| < |T/ 
for general /. Using ()8.13|) and (j8.14p . we have the bound 

(10.23) J (5"T"w + S^T^Huj)] ■ \S''T'^U\ ■ IS^T'IJl 

<(t}-3/2(||^3/2^ . (^'^T%)||ioc, + Wr^/^io ■ iS''T''Hu;)\\L^)E^AU] 

I/3|<1 

Pausing momentarily, we recall the definition of Af as in (|2.19bp and use Lemma 
EH ([STO]) . and the smallness oi eI(^[U] to derive 

(10.24) ^ l|r^'^T"+'^M^(i7)|| < ^ ||(r)|^''iT"i[/| ■ jS'^^T^^t/lll 

|/3|<1 |ail + |a2l<|o|+2 

ai + a2<a 

^ ||(r)S'"^T"i[7||Loo||5"2T"^[7|| 

|Ql| + |a2l<|o|+2 
oi+a2<o 
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Therefore, 

(10.25) J 7|w • (5°T"{> + ^"T"ijw)| • \S''T^U\ ■ |5'=T'^{/| 

Gathering estimates (|10.15|) . (ll0.16p . (|10.20p . ()10.2ip . and (|10.25|) . we have shown 
that each term in p0.14p is bounded by 



When we combine this result with (|10.13p . we have 



(10.26) dtE^,g[U] / '§"T"^;'5''T"Ai;' < {ty^'^ El[^g^^[U]E^,e+i[U]. 

a.<S 



Induction on 9 

We wiU complete the proof by proving estimates of the form 



(10.27) d,E,.jU] + - Yl \\yS'^T'^vr<{t)-'/'E,jU]El^'^^,[U] 

a+\a\<fj. 



and 



(10.28) dtE,,,+^[U] + - Y \\^S^r'^vr<{t)-'El/^yU]E,,,+,[U]. 

a<M + l 



via finite induction on 9 using (fTa9)) . (|lG.10p . and (|lG.26p . 
Applying (|10.26p with 9 = 1,2,...,^ and we have 



(10.29) dtE^^U] -vY. (Tf"«,Tf"Ai;) < {t)-^'^ El{l[U]E^,y[Ul 

|q|<M 

a+|c|<^ 



a+|a|</j 



The first line may be integrated by parts giving 

(10.30) d^E,AU\ + v Y l|VT"t;f <(t}-3/2iji/f[{/]iJ^^i[C/]. 

|a|<M 
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a+\a\<^L 



In the remaining lines we calculate by brute force using the commutation formula 
()3.10|) and integration by parts 

(10.31) V (5"'T"z;,S"^T"Az;) 



^+l°l<M 
o<9 



J2 / 5'''T"w'A5^T" 

-| = I<M i=0 ^-'^ 

J2 ||V5'^T"vf 



"+|c«l<M 



+|oi<M i=o ^-'^ 



V 



a+|a|<,, a+|c,|<,i 

a-1 



i E E(;)iiv^^t".| 



a+|a|<j, a+|Q|<M 



where 



(10.32) Ce = max Q 
Now we can write 

(10.33) + ^ ^ ||v^«T"i;f-^ ^ ||V5'^T%|| 

a+|c|<M a+|a|<M 

a<e a<e-i 

By a similar argument, for 1 < 6* < /i we use (|10.9I) to get 

(10.34) atK,e[C/] + ^ llV^'-T^f - ^ ^ ||V5''T%|| 

a+|c«|<« a+|Q|<K 
a<8 a<e-l 
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and (|10.10p to get 

(10.35) + ^ J2 \\VS-T"v\\' J2 l|V^"T"i;f 

a+\c.\<n a+\a\<n 
a<^l+l a<ti 

Because the dtE terms on the left of (|10.30p . (|10.33p - (|10.35p are not necessarily 
positive, we write each these estimates in its integrated form: 



(10.36) 



I 1^ Jo JQ 

(10.37) E,4U] + '^ f WVS'^T^vf-'^ rilV^"T««|| 

a<e a<B-l 

< C j\T)-^I^E'X,[U]E^.e+^[U]dT, 
Jo 

(10.38) E,,e[U]+'^ E TllV^'^T-^f ^ /* ||V^'^T%|| 

a+|Q|<K "'0 a+|c,|<K "'0 

a<e a<g-i 

<C f\Tr'El(^[U]E^,e+i[U]dT, 
Jo 

and 

(10.39) E^,,+,[U] + ^ E r ||V^'^T"i;f - ^ ^ f \\VS^T^ 



<G \\rr^E]ll{U\E^^^^mdT. 
Jo 

If we add (|10.37p with 6i = 1 to (|10.36p . we have 

(10.40) i?^,i[C/]+i?^,o[C/] + ^ E ri|V^''T«t,f + ^ ^ /*||VT"«|| 

^a+|c.|<,, -^0 |a|<p-^" 

a<l 



<c\ {T)'^l\E^ll\U\E^,^\U\+E^Jl\V\E,AV\)dT 
Jo 



which imphes 
(10.41) 



a+|Q|<,. 

3<1 



For = 2, 3, /i if we assume 

(10.42) i?^,,_i[C/] + ^ 5] t \\^~S''^''vf<\\r)-^l''E'll\lJ\E^^eW^ 



a+|c|<,. • 

3<e-i 
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then adding (|10.37p multiplied by OCe to (110. 42|) we have 
(10.43) E^,e[U]^eCgE^^e-i[U]+'^ ^ /* ||V^"T%||2 

< C f\T)-'/\El/'[U]E^,e+i[U] + eCgEl(^[U]E^Am)d^ 



which gives us 



(10.44) E,AU]+'^ J2 f\\yS^^''vr<f\r)-'/'El/^^,[U]E,,e+i[U]dT 



Therefore, by our inductive argument at = ^, we have 



(10.45) E,,,[U]+'^ Yl f \\^S^T-vf< f {rr'^'El/'^^,[U]E,,,[U]dT. 
A similar inductive argument for the high energy using (|10.38p shows 

(10.46) E^Am + 'i E f\\yS^r''vr<f\r)-'E,,o+i[U]El/'AU]dr 
^ „+i„i<„ Jo Jo 



:i+|c|<» 



for e = 0,1,2, At 6* = /i+ 1 we use (|10.39p to get 



\\S7S^T"v\\' < f\T)-^E,.,+,[U]El{l[U]dT. 
Jo 



(10.47) £;„,^+i[C/] + - 

a+|Q|<» 
a<,i+l 

Estimates (|10.45l) and (|10.47p are equivalent to estimates (|10.27p and (|10.28p which 
in turn imply global existence. □ 
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